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Abstract 

Azimuthal distributions of rj and tt" mesons emitted at midrapidity in collisions 
of 1.9 AGeV 5**Ni+5**Ni and 2 AGeV '"'Ca+""*Ca are studied as a function of the 
number of projectile-like spectator nucleons. The observed anisotropy corresponds to a 
negative elliptic flow signal for r; mesons, indicating a preferred emission perpendicular 
to the reaction plane. In contrast, only small azimuthal anisotropies are observed for 
tt" mesons. This may indicate that rj mesons freeze out earlier from the fire ball than 
pions. 

1 Introduction 

Heavy-ion collisions at incident energies of 1-2 AGeV (Bevalac/SIS energy regime) are a 
unique tool to study nuclear matter at high density and temperature. According to various 
theoretical model calculations pi, 0, 0] nuclear matter can be compressed to 2-3 times the 
normal nuclear density and heated to temperatures in the order of 100 MeV in this energy 
regime. In addition, a fraction of 10-30% of the participating nucleons is excited to short- 
lived resonance states, mainly A(1232) and A^*(1535) resonances, which subsequently decay 
via meson emission y, g. Thus the meson observables can provide information on the 
dynamical evolution of the resonance population in compressed and excited nuclear matter. 
The production and decay of the A(1232) resonance is responsible for the production of 
pions, which are the most abundantly produced secondary particles in the 1-2 AGeV energy 



regime P, Ri ^ ■ ^ clean signature for the excitation of a high lying resonance is the detection 
of 77 mesons M. They originate almost exclusively from the decay of the TV* (1535) resonance 
and thus are sensitive to the abundance of this resonance M, HQ] . After the initial production 
through resonance decays, pions and i] mesons strongly interact with the surrounding hot 
nuclear matter in the interaction region and in cold spectator matter. The interactions in 
hot nuclear matter are the result of a complex cyclic process of generation, absorption, and 
re-emission of mesons [0, M, while in cold nuclear matter absorption dominates. These 
processes will influence experimental observables like the meson abundance as well as the 
azimuthal angle distributions of mesons with respect to the reaction plane. 

Azimuthal anisotropy in the emission of pions in symmetric as well as in asymmetric 
non-central heavy-ion collisions is a clearly established effect and was observed in a wide 
range of beam energies from 0.1 A GeV to 160 A GeV ||l^ - |^. Similar to the anisotropy 
of baryons |gj| , the effect is usually discussed in terms of directed and elliptic flow |g^, Q . 
In a Fourier expansion of the azimuthal angle distribution of particles N{A(p) with respect 
to the reaction plane 






N{Aip) = vq \ 1 + 2 y Vn cos nAip 



n>l 



the first term (~ (cosAif)) corresponds to the directed flow, while the second term (^ 
{cos2Aip)) represents the elliptic flow pq|. At midrapidity the directed flow (cosAip) of 
emitted particles vanishes for symmetry reasons and only the elliptic flow is present p5| |. 

In the Bevalac/SIS energy regime the elliptic flow of participant baryons emerging from 
the collisions at midrapidity was found to be oriented perpendicular to the reaction plane 
(negative elliptic flow, {cos2Aip) < 0) |l^, |2^, ^ ^. This effect was interpreted as a 
dynamical squeeze-out of nuclear matter due to the build-up of pressure in the interaction 
zone between two colliding nuclei pO| . Similar to baryons, negative elliptic flow was observed 
for high transverse momentum neutral and charged pions emitted at midrapidity in 1 AGeV 
Au+Au collisions at SIS (GSI) [Q |l3- However, the observed pion anisotropy was not 
been attributed to the expansion of nuclear matter, but rather to the strong final state 
interactions of pions with cold spectator matter transiently concentrated in the reaction 
plane ("shadowing"). Consequently, the chance for pion absorption is higher in the reaction 
plane than out-of-plane |31 , 32 . The negative elliptic flow of high energy pions observed in 



the Bevalac/SIS energy regime indicates that they freeze-out while the spectators are still 
close to the participant zone Q . At ultrarelativistic energies the time required for spectator 
fragments to pass the reaction zone (passage time) is so short that the shadowing effect is 
reduced. In this case the geometry of the participant zone favors the preferential in-plane 
emission of pions p^, 34|, 35| and hence positive elliptic flow, which was observed recently at 



AGS |£1[ and SPS ||22, |23| en ergies. However already at SIS in-plane emission of low energy 
pions were observed, see p6J. Hence similar scenario can hold for low energy pions even at 
SIS energies. 

The elliptic flow of pions is well established and can qualitatively be described in the 
frame of transport model calculations p^, B2l p3, p5[ , but the elliptic flow of 77 mesons has 
not been observed before. While the absorption of pions in hot nuclear matter proceeds via 
the A(1232) resonance which decays dominantly by pion emission, only about 50% of all 
N*(1535) resonances excited by 77-meson absorption will reemit 77 mesons p(| . However, the 
mean free path of both mesons in cold nuclear matter is comparable 36, [37[. Therefore, a 
comparison of the 77- and 7r°-elliptic flow may yield information on the propagation of these 
mesons, as well as on the dynamics of the parent baryon resonances. 

Moreover, 77 and tt*^ mesons can be measured simultaneously in the same experiment via 
their two-photon decay channel (branching ratios: 39.3% and 98.8% for 77 and 7r°, respec- 
tively yi ) . Below we present results of the first experimental study of azimuthal distributions 
of 77 mesons emitted in collisions of 1.9 AGeV ssNi-fS^Ni and 2 AGeV ''•'Ca4-""*Ca nuclei 
and compare them with azimuthal distributions of vr*^ mesons in the same colliding systems. 




Figure 1: Schematic view of the experimental setup with TAPS and the KaoS 
Forward WaU. 



2 Experiment 

The experiments were performed at the Heavy-fon Synchrotron SIS at GSI Darmstadt. The 
experimental setup is sketched in Fig. |l|. In the first experiment a 1.9 AGeV ^^Ni beam 
with an intensity of 6.5 x 10^ particles per spill (spill duration 8 s and repetition rate 15 s) 
was incident on a ^*Ni target (502 mg/cm'^). In the second experiment, a "°*Ca target (320 
mg/cm^) was bombarded by a '*'^Ca beam with kinetic energy 2 AGeV and an intensity of 
5 X 10^ particles per spill (spill duration 10 s and repetition rate 14 s). 

Photon pairs from the neutral-meson decay were detected in the Two- Arm Photon Spec- 
trometer (TAPS) 1 38 . This detector system consisted of 384 BaF2 scintillators arranged in 



6 blocks of 64 modules with individual Charged Particle Veto detectors (CPV) in front of 
each module. The blocks were mounted in two towers positioned at 40° with respect to the 
beam direction at the distance of 150 cm. Three blocks were positioned in each tower at 
-1-21°, 0° and -21° with respect to the horizontal plane. In this setup, only neutral mesons 
around mid-rapidity j/cm were detected. The geometrical acceptance of TAPS for the tt*^ 
and 7] detection was roughly 1 x lO"'^. 

An in-beam plastic scintillator (BC-418) of 200 fim thickness was used to provide a time- 
zero signal for the time-of-flight measurements as well as to count beam particles. The plastic 
Forward Wall (FW) of the KaoS collaboration, see p3], comprising 380 plastic scintillators 
(BC408) was positioned 520 cm downstream of the target. The modules have a thickness 
of 2.54 cm with sizes 4x4 cm^ in the center, followed by 8x8 cm^ and 16x16 cm^ elements 
in the outer region. Four modules in the center of the wall are removed for a beam pipe 
made of a carbon fibre tube (diameter 70 mm, thickness 1 mm). The actual measurements 
were performed with 320 detectors to obtain a nearly azimuthally symmetric coverage. In 
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Figure 2: The phase space ( transverse momentum pt versus rapidity Y) for baryons 
from the '"'Ca+''°Ca at 2.0 A GeV coUisions calculated with FREESCO model @. The 
reaction detector covered the polar angles QLab from 14° to 30° (dotted lines). The solid 
lines indicate the region of the acceptance of the FW detector (6Lafc=0.7-10.5°). The 
rapidity is the laboratory rapidity. The center of mass rapidity is ycm=0.882. 



this arrangement the FW covered the polar angles from 0.7° to 10.5*^. Particles emitted 
in this angular range are predominantly projectile-like spectator nucleons, see Fig. 0. The 
FW provided the information on emission angle, multiplicity, charge and time-of-flight of 
protons and light charged fragments up to Z—8, see Fig. 0. 

The reaction detector, comprising 40 small plastic scintillators (NE102A), was positioned 
close to the target and covered the polar angles from 14° to 30°. Most of the particles 
emitted in this angular range are participant nucleons, see Fig. 0. For the selection of 
events according to the reaction centrality and for the estimation of the average number of 
projectile-like spectator nucleons the information from the reaction detector and the FW 
were used. The performance of these two detectors is shown in Fig. 0, where the total charge 
Zpw of the particles detected by the FW is plotted as a function of the charged-particle 
multiplicity M^eact in the reaction detector, measured in coincidence with two neutral hits 
in different blocks of TAPS. 

Events were selected according to the reaction centrality by requiring equal statistics of 
mesons in each bin, ranging from peripheral reactions with low multiplicity Mreact to central 
reactions with high multiplicity Mreact ■ The total charge Zpw of particles detected by the 
FW allowed us to estimate the mean number of projectile-like spectators (Agp) for each 
studied bin in Mreact as determined by the reaction detector. We used the relation 



i^sp) — {ZFw}Aproj/Z^^ 



proj j ^proj t 



(2) 



where Ar^m-i and Z^^-o^ are the mass number and charge of the projectile, respectively. The 
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distributions of the total charge Zpy/ are shown in Fig. ^ for both the Ni+Ni and Ca+Ca 
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Figure 3: The dependence of energy loss dE/dx as a function of particle velocity /3 is shown 
for fragments detected in a single FW module. The data are from the ^^Ni+^^Ni reaction 
at 1.9 AGeV. 
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Figure 4: The total charge T^fw of particles detected in the FW as a function of the charged- 
particle multiplicity Mreact, measured in coincidence with two neutral hits with an energy above 
90 McV in two different TAPS blocks. The data are from the ^^Ni+^^Ni reaction at 1.9 AGeV. 
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Figure 5: The distribution of the total charge of particles detected in the FW for different bins 
in charged-particle multiplicity Mreact giving the reaction centrality. (a) - (c) for the experiment 
^^Ni+^^Ni at 1.9 AGeV; (d) - (e) for the experiment ''"Ca+''°Ca at 2 AGeV. The arrows indicate 
the mean values of the total charge. 



Reaction 


b»Ni+b8Ni at 1.9 A GeV 


^Ca+natQg^ at 2 A GeV 


-LVlreaci 


2-6 


7- 10 


> 11 


1-3 


> 4 


(Asp) 


51 


44 


37 


34 


28 



Table 1: The mean number of projectile-like spectators (Asp) for studied bins in the 
charged-particle multiplicity Mreact measured by the reaction detector. 



collisions. The resulting values of the mean number of projectile-like spectators {Agp) are 
listed in the Table IT]. The systematic error of the values (Agp) was found |El| to be less than 
4 units. 

Neutral-meson candidate events were selected by special triggers exploiting the kinemat- 
ical constraints of meson decay into photons. For each TAPS block the hit information for 
the BaF2 signals above a certain threshold energy and for the corresponding CPV modules 
were fed to a multiplicity/pattern unit. This unit performed an on-line selection of events 
induced either by a charged particle (if the BaF2 scintillator and the corresponding CPV 
module have fired simultaneously) or by a neutral particle (if the BaF2 scintillator has fired 
without a corresponding CPV signal). For the analysis of the present data the following 
triggers were chosen: 
For the reaction ''*'Ca-(-""*Ca at 2.0 A GeV the event-selection trigger required 

a) hits with an energy E>15 MeV in any two TAPS blocks for the 7r° measurement, and 

b) at least one neutral hit in each tower with an energy E>90 MeV for rj mesons. 
For the reaction ^*Ni-|-^^Ni at 1.9 A GeV the trigger required 

a) one neutral hit with an energy E>90 MeV and one hit with an energy E>15 MeV in any 
two TAPS blocks for a tt" candidate, and 

b) two neutral hits with an energy E>90 MeV in any two TAPS blocks for selection of rj 
meson candidate events. 



In addition, all triggers mentioned above required the coincidence with a signal from the 
reaction detector and the FW. 

In both reactions downscaled data were taken corresponding to minimum bias trigger re- 
quiring coincidence between the reaction detector and the FW. 

3 Data Analysis and Simulation 

3.1 Reaction-plane determination 

The study of event anisotropy requires the determination of the orientation of the reaction 
plane. With respect to this plane the preferred emission of some subset of particles can be 
analysed. The reaction plane was determined by a modified version p2[ of the transverse- 
momentum method [Esj . For each event the reaction plane was defined by the incident beam 
direction and the vector Q, which is the weighted sum of the transverse momentum vectors 
of all charged particles detected in an event. Since almost all particles detected by the FW 
have nearly the same velocity (see Fig. |[), the transverse velocity vector v-j^ will be oriented 
parallel to the position vector r^ of the particle k in the x-y plane (perpendicular to the 
beam direction) of the FW. Here it is assumed that the beam position in the plane of the 
FW is (0,0). Therefore, in our analysis the vector Q is defined as 



Q-E 



M 

^fcr^> (3) 



where the sum runs over all M particles detected by the FW in the event, r^ is the position 
vector of particle k in the x-y plane and uJk is the weight factor, which depends on the 
rapidity of the emitted particle. Since the coverage of the FW excludes nearly all particles 
with j/fc < Hem the position vectors in Eq. 3 will be summed with a positive weight Uk. In 
order to provide the best reaction plane resolution the weight factor ujk was chosen to be 
ujk=Zk, where Zk is the charge of particle k detected by the FW. 

Since the coverage of the FW does not overlap with the TAPS spectrometer (used for 
the neutral meson reconstruction), the analysis of azimuthal angle distributions of 77 and tt" 
mesons with respect to the reaction plane is free from autocorrelation effects. The analysis 
of azimuthal anisotropics is rather sensitive to different kinds of experimental biases, which 
could simulate an event anisotropy. Such artificially created anisotropics can be removed 
by requiring the distribution of the reaction plane to be isotropic for unbiased events. The 
off-line data analysis revealed that the actual beam position in the plane of the FW with 
respect to the geometrical center of the FW varied with time. Consequently, the distribution 
of the reconstructed reaction-plane angle is not flat. This effect required a correction of the 
vector f*"'"'' = fk — Toff for each particle in an event. The offset values {{xoff), iUoff)) 
were determined by averaging over event samples in fixed time intervals. The required 
corrections {{xoff), {Uoff)) stayed below 1 cm. As an example of this procedure, Fig. ra 
shows the distributions of the reaction-plane angle ($_r) for peripheral (2 < Mreact < 6) 
Ni-|-Ni reactions at 1.9 AGeV before and after including the beam-offset correction. The 
corrected distribution (Fig. m) is free from significant distortions which could influence the 
study of azimuthal anisotropics. 

3.2 Reaction-plane resolution 

Because of finite multiplicity fluctuations, the azimuthal angle ^r of the vector Q can 
differ from the azimuthal angle of the true reaction plane (j^true by a deviation A(j)pi — 
(f'true — ^R- For Corresponding coefficients vf"^^ and v^^'^^ of Fourier expansion of azimuthal 
angle distribution of particles with respect to the "true" reaction plane then holds 1;^'""'^ = 
vi/{cosA(j)pi) and vl™"^— f2/(cos2A(/)p/). 
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Figure 6: The distributions of the reaction plane angle ^r for peripheral ( 2< Mreact <6 
) ^*Ni+^*Ni reactions before (a) and after (b) the beam-offset correction. 




Figure 7: The variation of {cos nA(/)S"'*) with the parameter x, calculated from Eq. 5 for 
the first two harmonics n=l and n=2. 
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Figure 8: The distribution of the relative azimuthal angle A$i,2 between the reaction 
planes Qi and Q2 of two subevents. The right part of the figure shows data for a well 
defined reaction plane with | Q |> 3. The data are from peripheral (2< Mreact <6) 
^^Ni+^^Ni reactions. 



For the determination of these corrections we follow the procedure as detailed in ref. 
|44| , [45[ . The distribution of A(/)p; can be presented as a function of a single dimensionless 
parameter x which measures the accuracy of the reaction-plane determination. For a large 
sample of events within the same centrality window, according to the central limit theorem, 
the fluctuations of the length Q =| Q | of the reaction- plane vector around its average value, 
(Q), are distributed as a Gaussian: 



dN 



1 



QdQdAippi TTcr- 



:exp 



\Q-{QW 
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jexp 



Q2 + (Q)2 _ 2Q{Q)c 



(4) 



The fluctuations with standard deviation a are assumed to be isotropic since the azimuthal 
anisotropies are small. The dimensionless parameter x={Q)/o- scales with the particle multi- 
plicity M like y/M. Integration of Eq. 4 over both A0p/ and Q yields the value {cos nA<j)pi), 
which will be used later for the correction of the measured azimuthal anisotropies. According 
to the "OUitrault method" we find: 



{cos nAif)pi) 



{cos nA<j,^l'^') 
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-XV2 



/.- 



(5) 



where Ik is the modified Bessel function of the order k. 

Fig. shows the variations of the first two correction coefficients for n=l (correction for 
the directed fiow signal) and n=2 (correction for the elliptic fiow signal) with x- 

In order to determine the parameter x ^^ randomly divided the hits in each event into 
two subgroups containing each one half of the number of particles. For each of the two 
subevents one can construct according to Eq. 3 the two independent vectors Qi and Q2, 
respectively, and extract the angle A$i_2='I'it',,i — ^r,2 between the two vectors. Fig. || 
shows the distribution of the relative azimuthal angle A<I>i 2- 

The ratio of events with | A$i_2 |> 90° to the total number of events allows to determine 
the parameter x from Eq. 6: 



iV(90° < A$i,2 < 180") _ exp{-x^/2) 
iV(0O < A$i,2 < ISO") ~ 2 



(6) 



Reaction 


b»Ni+^«Ni at 1.9 A GeV 


4uca+"'^*Ca at 2A GcV 


Mreact 


2-6 


7- 10 


> 11 


1-3 


> 4 


(^pl 


43° 


46° 


50° 


52° 


55° 


{cos2^4>T') 


0.37 


0.32 


0.26 


0.22 


0.19 


{cos2A4,X'') 


0.33 


0.28 


0.23 


0.18 


0.16 



Table 2: The width (standard deviation api) of the relative azimuthal-angle distributions, 
the reaction plane resolution for the elliptic-flow signal (n=2) determined by "OUitrault 
method" {{cos2/^(jPl^'^^)) and by the adapted transverse-momentum method ( (cos2A(/)5"")) 
(see text) for studied bins in hit multiplicity Mreact of the reaction detector. 



Alternatively, these corrections can be estimated by the adapted |Q - |Q transverse- 
momentum method of Danielewicz and Odyniec p3| . First, we determined the average of 
the weighted unit vector of transverse velocity ly^^""^) of each charged particle, projected 
on the reaction plane. The orientation of the reaction plane is estimated by the vector Q™, 
determined here by excluding the considered charged particle in order to remove autocorre- 
lations: 



{V^ 



QT 



Ql 



QT^ 



M 



(7) 



Then, we calculated the average of the unit vector of transverse velocity {V^'^^'^) projected 
on the "true" reaction plane. 
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(8) 



The corrections due to finite reaction-plane resolution for the directed-flow signal (n=l) and 
the elliptic-flow signal (n=2) can be determined from the following equations: 



{cosA(f>: 
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{V^' 
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(cos2A(/i, 
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= 2 . (cos'AcP^r) 



- 1 



(9) 



Our further analysis was restricted to a sufficient vector length | Q |> 3 in order to reject 
the most central events lacking spectator flow. This selection rejects 25% of all registered 
events. The resulting values of the correction coefficient for the elliptic-flow signal due to 
finite reaction-plane resolution (cos2A0p;), determined by both methods, are given in Table 
for each bin in hit multiplicity Mreact of the reaction detector for both systems studied. 
The systematic error of these values was estimated from the difference |(cos2A(/)S''") — 
{cos2A(j)E°'^)\ and is found to be less than 15%. We also studied the influence of the cut 

I Q |> 3 on these values. The removal of the cut leads to decrease in the correction value, 
see left panel of Fig. g. and Eq.6. But simultaneously also deduced values of V2 decrease 
and hence the Wj'^"'^ values remain unchanged. 

For comparison with previous studies we calculated also reaction plane resolution api — 
cri.2/2, where standard deviation ai^2 is extracted from a Gaussian fit to the dN/dA^i^2 
distributions. Reaction plane resolution varies between 43° and 55° depending on the re- 
action centrality and the colliding system (see Table g) . These values agree with published 
data from studies of charged-baryon flow in similar colliding systems [^3[ ^. Moreover, 
the the directed flow of charged baryons in the target-like rapidity region was observed for 
both colliding systems (see Appendix A). This verification demonstrates the quality of the 
reaction-plane determination. 
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3.3 Photon-particle discrimination 

EfScient photon-particle discrimination in TAPS can be achieved by the combination of 
three methods |l| - ||: 

1) The anticoincidence between the BaF2 and the corresponding CPV was required to 

discriminate charged particle hits. 

2) The timc-of- flight analysis allowed to separate photons from massive particles (with non- 

relativistic velocities) . 

3) The pulse-shape of the BaF2 energy signal was analyzed. The scintillation light in the 

BaF2 crystal exhibits two components with different time constants r, one fast (t=0.6 
ns) and one slow (t=620 ns). The intensity ratio of fast to slow component depends 
on the ionization density of the incident radiation: it is larger for leptons (electrons 
following photon interactions) than for other charged particles (p,d,t,a...). The pulse- 
shape (PSA) parameter is defined as the angle ^psA — s^'^^ctan^Eshort/Etotai), where 
^short is the fraction of energy obtained by integrating the analog signal within the 
first 30 ns, and Etotai is the total energy obtained by integrating the analog signal 
within 2/is. Fig. ya presents a two-dimensional graph of Eghort versus Etotai for all 
BaF2 modules, and the distribution of the corresponding PSA parameter is shown in 
Fig. ^. The separation between photons and charged particles is obvious. 



3.4 Shower reconstruction by cluster analysis 

Hadrons, like protons and charged mesons, lose their energy by excitation and ionization. 
On the other hand, high energy photons impinging on TAPS crystals generate an electro- 
magnetic shower which fires several neighbouring detector modules. Therefore, the recon- 
struction algorithm to separate photons and hadrons and to recover position and energy 
of the incident photon is based on defining a shower as a continuous cluster of individual 
responding detectors (52| |5J, |5^, Q. In the first step, the fired detector modules are sorted 
according to the deposited energy. The combination to clusters starts with the detector 
which indicates the maximum energy deposit in a block of crystals above a certain en- 
ergy threshold Fihigh (typically ~20 MeV). This so-called central detector is combined with 
all directly neighbouring detectors (maximum 6 modules in the first ring of neighbours). 
However, these modules are only accepted in the cluster if the time difference between the 
neighbouring and the central module is within the allowed time window and an energy depo- 
sition above the lower threshold Eiow (typically ~3 MeV) is encountered. If the sum energy 
of this minimal cluster exceeds a threshold of typically 400-500 MeV or if the fraction of 
the energy deposited in the central module is lower than 60-70 % of the sum energy, then 
the second ring of neighbouring detectors (up to 12 modules) is inspected for addition to 
the existing cluster. All allocated detectors are then removed from the detector hit pattern 
and the cluster search continues with the remaining highest-energy module. This procedure 
limits the cluster size to a maximum of 19 detector modules which is, however, sufficient to 
reconstruct even the highest energy photons encountered in this experiment. 

In the second step the characteristics of all modules within a cluster were compared to 
distinguish different types of showers. The cluster is called a photon- like cluster if: 

1. the hit in the central detector is neutral, i.e. no CPV signal was found in the correspond- 

ing CPV module or the CPV's of the neighbouring detectors. 

2. the pulse shape corresponds to the hit of a photon (see Fig. 0). 

3. the time-of-flight is within the prompt time window. The neighbouring energy deposition 

is added only if the neighbours also have photon-like characteristics. 
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Figure 9: Pulse-shape analysis of BaF2 signals: a) graph oiEghort versus Etotai for the 7r° trigger 
in Ni+Ni reactions at 1.9 A GeV; b) distribution of the pulse-shape parameter ^psA for all 384 
BaF2 modules. In this case the pulse-shape analysis was performed only for Etotai >30 MeV. 
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The average particle occupancy in TAPS reaches values of 7% per module. This leads to 
a finite probability that the clusters in the TAPS block are connected or partly overlapping. 
However, in general connected or overlapping clusters have a second maximum and there 
will be more energy deposited in the second ring as is expected from a single hit. The 
clusters with a second maximum were removed from the further analysis. The energy of the 
cluster E results from the sum of energies Ei deposited in the individual detectors, and the 
direction of the incident photon r is reconstructed from an energy- weighted sum over the 
vectors r^ pointing from the target to each responding detector from the cluster [p7|. 



E^y^E, and r = ^^'^' ^jhere W^ = MAX ( 0, 



Ei' 
Wo + In-^ 



(10) 



Parameter Wq controls the importance of the positions of the neighbouring modules with 
respect to the position of the center module. 

3.5 Neutral-meson reconstruction 

Neutral pious and ?7-mesons can be detected by the measurement of their two-photon decay, 
which still carries the signature of the meson in its invariant mass. For each pair of detected 
photons in a given event we calculated the invariant mass Mpair and the momenta Ppair 
using the following relations: 

*^pair = 2^1^2(1 -COsOia) and Ppair^Pl+P2, (11) 

where Ei, pi and E2, P2 are the energies and momenta of the corresponding photons, and 
012 is the opening angle of the photon pair. We analysed only neutral mesons in a narrow 
rapidity window (|y — ycm\ < 0.1). The resulting invariant-mass distributions (Fig. |o|- ^ 
reveal peaks centered at values of the 7r°- and yy-meson rest mass. These peaks are located 
on top of a large combinatorial background which originates from uncorrelated photon pairs. 
The high particle multiplicity in heavy-ion collisions, limited acceptance and the detector 
efficiency usually lead to very low signal-to-background ratios, especially at lower transverse 
momenta. Therefore, the precise knowledge of the shape and magnitude of the combinatorial 
background is required for a reliable evaluation of the neutral-meson peak contents. The 
analysis quality can be improved by increasing the purity of the photon sample, i.e., by 
reducing the combinatorial background. Therefore, all cuts in the pulse-shape analysis and 
cluster reconstruction scheme were optimized for a maximum signal-to-background ratio 
in the invariant-mass spectrum. The resulting signal-to-background ratios for the studied 
bins in reaction centrality M^eact and transverse momentum pt are listed in the Table ^ 
for both reactions studied. If the background has a monotonous falling or rising shape it 
will not be difficult to approximate the shape of the background with an analytic function. 
However, usually the artificial structure in the invariant-mass spectra created by limited 
detector acceptance could not be reproduced by any analytical fit of the background. 

3.6 Method of event mixing 

Since the combinatorial background consists of pairs of uncorrelated particles, the natural 
way of its reconstruction is the use of the so-called " event-mixing" technique by combining 
photons randomly selected from different events [gSJ - |Q. This method is based on the 
assumption that photons from different events are not correlated. Therefore, the mass signals 
( tt", 77-mesons) are absent when we construct a two-photon invariant mass distribution by 
taking photons from different events. This procedure automatically takes into account the 
detector acceptance and the detection efficiency. Such combinations can be formed in large 
numbers, so that the combinatorial background can be determined with high statistical 
accuracy. Then the background distribution is normalized outside of the meson-mass window 
to the measured invariant-mass distribution and subtracted from the latter. The meson-peak 
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Figure 10: Results of the invariant-mass analysis of photon pairs. The upper frame 
shows the invariant-mass spectrum which corresponds to the 77 trigger in the experiment 
^^Ni-|-^*Ni at 1.9 AGeV. The combinatorial background (dotted line) was determined by 
event mixing. The lower frame shows the invariant-mass distribution after background 
subtraction and demonstrates the quality of the background determination. 

contents can be extracted from the background-free invariant-mass spectra by integrating 
within the meson-mass window. 

However, for a perfect description of the combinatorial background special care has to 
be taken to ensure that the data and background distributions are generated for the same 
class of events with very similar event characteristics. Therefore, only photons from the 
same event class may be mixed. An event class is defined by the analysed trigger TRG and 
the photon multiplicity M-^ found in the event. Including the constraint on M-,, in the event 
classification is important as the signal-to-background ratio (S/B) depends on M^ |6l]. The 
phase-space distribution of the mixed photon pairs must correspond to the measured ones. 
The phase space covered by TAPS allows the identification of neutral mesons emitted in 
a narrow window around midrapidity only. As a result, the transverse momentum of the 
photon pair can classify the phase space allocation. Therefore, if one divides the photon pairs 
in the classes according to TRG and M^ into groups according to their transverse momenta 
Pf, an identical phase space allocation for mixed and measured pairs can be expected. As 
an example, Fig. O shows the invariant-mass distributions for different P(-classes and the 
photon-multiplicity class M.y=2 (tt^ trigger for the experiment Ni-fNi at 1.9 A GeV). The 
normalized combinatorial background is shown as a dashed line together with the histogram 
for the experimental data in each case. The spectral form of the distributions varies strongly 
with the transverse momentum which refiects the opening angle acceptance of TAPS. 

Subsequently, the mixed photon pairs were selected to have an opening-angle distribu- 
tion close to the opening angle of the original pair. This selection can be implemented by 
an additional classification according to the TAPS block combination. However, the initial 
concept of a good reproduction of the phase-space occupancy of the invariant-mass distri- 
bution in the combinatorial background by the method of event mixing can still suffer in 
some cases. Besides the true correlations due to the neutral-meson decay there are some 
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Figure 11: Results of the invariant-mass analysis of photon pairs. The upper frame shows the 
invariant-mass spectrum which corresponds to the 77 trigger in the experiment Ca-|-"'^*Ca at 
2.0 AGeV. The combinatorial background (dotted line) was determined by event mixing. The 
lower frame shows the invariant-mass distribution after background subtraction. 
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Figure 12: Results of the invariant-mass analysis of photon pairs. The upper frame shows the 
invariant-mass spectrum which corresponds to the 7r° trigger in the experiment '*°Ca-|-"°*Ca 
at 2.0 AGeV. The combinatorial background (dotted line) was determined by event mixing. 
The lower frame shows the invariant-mass distribution after background subtraction. 
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Reaction 

^react 


i^«Ni-F^«Ni at 1.9 A GeV 
2-6 7-10 > 11 


4'^Ca+ 
1- 3 


'"^*Ca at 2 A GeV 
> 4 


Pt (MeV/c) 
77 0-600 


3.9% 


3.1% 


2.5% 


8.9% 


4.5% 


tt" - 200 
200-400 
400-600 
600-800 


9.5% 
21% 
61% 
83% 


6.2% 
13% 
33% 

47% 


4.3% 
9.2% 
24% 
36% 


15% 
36% 
99% 
140% 


10% 
24% 
65% 
110% 



Table 3: The signal-to-background ratio (S/B) for the tt" and 77 reconstruction for different 
bins in multiplicity M^eact and transverse momentum p^. 
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Figure 13: Invariant-mass spectra corresponding to the tt" trigger in the experiment 
^®Ni-|-^*Ni at 1.9 AGeV for different bins in transverse momentum pt (range in MeV/c is 
noted in each spectrum) and photon multiplicity M^=2. The combinatorial background 
was determined by event mixing and is shown as a dashed line together with the histogram 
for the experimental data in each case. 
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Figure 14: Azimuthal-angle distributions of inclusive photons detected by TAPS from pe- 
ripheral (2< Mreact <6) Ni+M reactions at 1.9 AGeV, averaged over photon transverse 
momentum. The solid line is a result of the fit by the first two terms of the Fourier expansion 



additional correlations due to reaction dynamics or detector effects. These correlations are 
also destroyed by event mixing and therefore are not removed from the data. 

In a mixed event two photons can be very close to each other in space. In a real event, if 
two photons are too close to each other, they are not resolved, and are counted as one photon. 
This effect was reduced by storing photons only when they belong to different detector 
blocks as in the experiment. This leads to a small under-estimation of the background from 
mixing events, but only for very small invariant masses (below 50 MeV). The photon pairs 
from meson-peak regions (7r*',?7) will be more abundant than photon pairs from the regions 
outside of it |pq]. Therefore, if the combinatorial background will be generated from all 
photon pairs from the meson-peak regions, we will overestimate the contribution from these 
photons. The ratio of the yield of the meson peak to the yield of the total invariant mass 
S/T=Signal/ (Signal -|-Background) was used to suppress photon pairs from the resonance 
regions for each event class (TRG, M-y, pf). 

We observed weak azimuthal anisotropy found in the yield of inclusive photons detected 
by TAPS. As an example, Fig. |lj shows the azimuthal angle distributions of inclusive pho- 
tons detected by TAPS with respect to the reaction plane from peripheral (2< Mreact <6) 
Ni-|-Ni reactions at 1.9 A GeV averaged over photon transverse momentum. The Fourier 
analysis of this distribution shows the presence of weak positive elliptic flow of inclusive 
photons V2 — .01, indicating a preferred emission of photons in the reaction plane. Further- 
more, the detailed analysis revealed that the magnitude of the combinatorial background 
is also weakly dependent (^2^ — .01) on the azimuthal angle At/? = (l>pair — ^r of meson 
emission relative to the reaction plane. As the main source of photons in our case is due to 
tt" decay, both these anisotropics can be attributed to the anisotropy in the pion emission. 
As the integrated pion yields are dominated by pions with low transverse momentum pt, 
the elliptic flow of low pt neutral pions is also expected to be weak and positive. The per- 
formed Monte-Carlo simulation confirmed this assumption. The simulation also shows that 
the correction for this effect can be done by normalizing the combinatorial background for 
each bin in azimuthal angle A(/3 separately. 
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Figure 15: The extraction of the ry-meson peak contents in the interval 0<pt < 400 MeV/c in 
the 12 equal intervals in azimuthal angle Ai^j = (l)yair — ^r of meson emission relative to the 
reaction plane ^r. The fitted line is an asymmetric Gaussian distribution [p3|. The number 
of 77 mesons for each bin in azimuthal angle </> was obtained by integrating the background- 
free invariant-mass distribution within the meson-mass window: 470 MeV/c^ <m^^ < 570 
MeV/c2. The reaction is 40Ca+'"'*Ca at 2 AGeV. 
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Figure 16: Azimuthal-angle distribution of 77 mesons emitted in '"'Ca-|-""*Ca collisions at 2.0 
AGeV averaged over multiplicity Mreact and transverse momentum pf The solid line is the 
result of a fit by the first two terms of the Fourier expansion. 
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In our analysis we used 12 equal intervals in azimuthal angle A(^ = (j^pair ^ ^_r of meson 
emission relative to the reaction plane '^r. The normalization factor was found to be a 
smooth function of A(/9. 

In order to extract the azimuthal-angle distributions of mesons the following algorithm 
was used: In a first step, the calculated invariant-mass distributions for each studied trigger 
and reaction-centrality window were divided into several sections according to the photon 
multiplicity M^, transverse momenta pt, TAPS block combination and azimuthal angle A(/3. 
Then the combinatorial background was generated by the outlined scheme of event mixing, 
normalized outside of the meson-mass windows and subtracted from the invariant-mass 
spectra for each section separately. Then the background-free invariant-mass spectra from 
each section were summed. As an example. Fig. y_5| shows the resulting background-free 
invariant mass spectra in the 77 meson mass region for 12 equal intervals in Aiy9 summed over 
the reaction centrality M^eact and transverse momentum p^ . The 77 meson peak content was 
obtained by integrating over the mass window 470 MeV/c^ <M^^ <570 MeV/c^, and for 
the tt" mesons we used the mass window 100 MeV/c^ <M^-y <150 MeV/c^. The resulting 
azimuthal-angle distribution of 77 mesons is shown in Fig. |l^. 

The uncertainty in the meson-yield extraction is determined by the statistical errors in 
the contents of the total mass distribution in between the integration limits and the number 
of entries in the (scaled) combinatorial background distribution. The combined statistical 
error is taken to be the square root of the quadratic sum of these contributions. 



3.7 Simulation 

As the main source of combinatorial background in our case are uncorrelated photon pairs 
which mainly originate from i:^ decay, the weak anisotropy of the background can be at- 
tributed to the anisotropy in pion emission. To test this assumption and to check the accu- 
racy of the proposed procedure of event mixing a Monte-Carlo simulation was performed in 
which the emission of tt" and 77 mesons is governed by a thermal distribution 

dN 



, ex const ■ X ■ y/x'^ - 1 • Ee^^'^/^'^ (12) 

dE 

where E and x = E/M arc the total energy of the meson ( M is the meson mass ) and its 
reduced total energy in the center-of-mass system, respectively. The meson ( tt", 77 ) source 
parameters, i.e. temperature T and meson production probability per participant were 
chosen consistently with those measured ||6J]. The polar-angle distribution of the emitted 
mesons is isotropic in the center-of-mass system. 

The azimuthal-angle distribution iV(A(p) of the emitted mesons with respect to the 
reaction plane is not isotropic and is given by : 

7V(A(/j) = const • (1 + 2 • 7;f "'^ • cos(2A(^)), (13) 

where the parameter ti|'""'^ reflects the strength of elliptic flow. The anisotropy in the emis- 
sion of rj mesons was assumed to be 7;2''"'^=-0.2, close to the experimental value. The 
simulations were performed with several values of the anisotropy in the emission of pions 
(—0.12 < v^^^ < -1-0.12), see below range of values actually observed in the experiment. Ge- 
ometrical acceptance, on-line trigger schemes and realistic detector resolution effects (time- 
of-flight and energy resolution of BaF2 modules) were taken into account. We used the 
identical procedure of analysis for the simulated data as for the experimental ones. The 
statistics in the simulated data was chosen close to that of the experimental data. The 
results of the simulation show that 

• the azimuthal anisotropy of the combinatorial background is an indication of the col- 
lective flow of pions, and its sign and magnitude are determined by the anisotropy of 
the pionic flow. Isotropy in the emission of pions leads to an azimuthally isotropic 
combinatorial background. 
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Figure 17: Azimuthal-angle distributions from a Monte-Carlo simulation, a) for tt*' from a pion 
source with strong positive elliptic flow; b) for 77 mesons from the 77 source; c) for r/ mesons 
reconstructed by the proposed scheme of event mixing. 
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Figure 18: Azimuthal angle distributions from a Monte-Carlo simulation, a) for 7r° from a pion 
source with a strong negative elliptic flow; b) for 77 mesons from the rj source; c) for rj mesons 
reconstructed by the proposed scheme of event mixing. 
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Figure 19: Azimuthal-angle distribution of 7r° mesons emitted in ^^Ni+^^Ni reactions 
at 1.9 AGeV. The panels correspond to different bins in reaction centrality M^eact and 
transverse momenta p* . 



• the extracted values for the elliptic-flow signal of 77 and 7r° mesons are, within statistical 
errors, comparable with source parameters; the proposed method of event mixing works 
correctly in case of different anisotropics of the pion source. As an example we present 
the results of such simulations for a pion source with strong positive (see Fig. 17) and 
negative elliptic flow (see Fig. |8|) 

Based on these simulation we estimate a systematical error in the parameter V2 for the 
experimental data, which arc influenced by much smaller mean pion anisotropy and hence 



also combinatorial background anisotropy v. 



bg 



-0.01, to be below 2%. 



4 Results 



We present in Fig. O - 22 the resulting azimuthal yields of 77 and tt" mesons for different 
bins in Mreact and transverse momentum pt for both systems studied. Because of low 
statistics for rj mesons we present the pt dependence (see Fig. 20) without a selection on 
the reaction centrality Mreact, and the reaction centrality dependence without a selection 
on the transverse momentum p^ (see Fig. 21). We fitted the azimuthal yields of rj and tt*^ 
mesons by the first two terms of a Fourier expansion in the azimuthal angle: 



iV(A^) 



2tt 



( 1 -f 2vicos{Aip) + 2v2Cos{2Aip) ) . 



(14) 



The first coefficient (vi ) is used to parametrize the in-plane emission of the particles parallel 
(vi > 0) or antiparallel (vi < 0) to the impact-parameter vector (directed flow), whereas the 
second coefficient (W2) quantifies an elliptic flow, which is negative for out-of-plane emission 
and positive for an in-plane emission [E5l u& . The extracted values of vi are zero within the 
error bars, as should be expected since we study symmetric colliding systems at midrapidity, 
see Tables B-pL The resulting values of the parameter V2 for 7/ and tt'^ mesons are given in 
Table | and § 

The obtained V2 coefficients were corrected for the reaction-plane resolution as V2^"'^ = 
V2/{cos2A(f)pi), where Acjjpi — cjjtrue — ^r is the fluctuation of the azimuthal angle of the 
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Figure 20: Azimuthal angle distribution of 7r° mesons emitted in ^*'Ca-|-'*''Ca reactions 
at 2.0 AGeV. The panels correspond to different bins in reaction centrality Mreact and 
transverse momenta pt . 



Reaction 
Pt (MeV/c) 


0-200 


h^'^Ni at 1.9 i' 
200-400 


tOeV 
400-600 


*"Ca+"''*C 
0-200 


a at 2 AGcV 
200-400 


Vl 


-0.08±0.07 


-0.03±0.05 


-0.07±0.08 


0.01±0.04 


0.02±0.04 


V2 


0.10±0.06 


0.11±0.05 


0.10±0.05 


0.11±0.04 


0.10±0.04 



Table 4: Parameters vi and V2 (not corrected for the reaction-plane resolution) for 77 mesons 
deduced from the experimental azimuthal angle distributions (see Fig. UV) for several inter- 
vals in transverse momentum pi , without a selection on the reaction centrality. 



Reaction 

^react 


^«Ni-|-^«Ni at 1.9 AGeV 
2-6 7-10 > 11 


«'Ca+""*Ca at 2 AGeV 
1-3 > 4 


Pt (MeV/c) 


Vl for rj mesons 


0-600 


0.01±0.06 


0.06±0.10 


0.11±0.16n 


0.05±0.06 


-0.02±0.08 




Vl for n'^ mesons 


0-200 
200-400 
400-600 
600-800 


-0.021±0.026 
-0.016±0.017 
-0.017±0.032 
-0.012±0.021 


0.011±0.022 
0.016±0.017 
0.010±0.022 
0.016±0.017 


-0.031±0.042 
-0.011±0.021 
0.006±0.017 
0.016±0.032 


-0.021±0.022 
0.006±0.011 
-0.005±0.011 
0.006±0.012 


0.021±0.022 
-0.011±0.012 
-0.005±0.011 
-0.016±0.021 



Table 5: Parameters Vi (not corrected for the reaction-plane resolution) for tt*' and 77 mesons 
deduced from the experimental azimuthal distributions for several intervals in reaction cen- 
trality. 
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Figure 21: Azimuthal-angle distribution of 77 mesons with respect to the reaction plane for 
different bins in transverse momentum pt and averaged over muhiphcity Mreact'- (a) - (c) 
for the experiment ^^Ni+^^Ni at 1.9 AGeV and (d) - (e) for the experiment '*°Ca-K"''*Ca at 
2 AGeV. 



Reaction 

^react 


^«Ni 
2-6 


+^»Ni at 1.9 A 
7- 10 


GeV 

> 11 


'"'Ca+""*C8 
1-3 


I at 2 AGeV 
> 4 


(Asp) 


51 


44 


37 


34 


28 


Pt (MeV/c) 


V2 for 77 mesons 


0-600 


-0.08±0.04 


-0.12±0.05 


-0.10±0.06 


-0.09±0.03 


-0.09±0.04 




V2 for tt" mesons 


0-200 
200-400 
400-600 
600-800 


0.033±0.017 
0.005±0.007 
-0.004±0.005 
-0.031±0.012 


0.021±0.019 
0.007±0.009 
0.008±0.011 
0.012±0.013 


0.026±0.031 
0.011±0.016 
-0.018±0.016 
0.022±0.021 


0.012±0.016 
-0.011±0.006 
0.006±0.007 
0.005±0.012 


-0.006±0.021 
-0.005±0.011 
0.016±0.007 
0.007±0.012 



Table 6: Parameters V2 (not corrected for the reaction-plane resolution) for 7r° and r] mesons 
deduced from the experimental azimuthal distributions for several intervals (a-e) in Mreact- 
The corresponding mean number of projectile-like spectators (Asp) are indicated (see text). 
For TT*^ mesons the selected intervals in transverse momentum pt are indicated as well. 
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Figure 22: Azimuthal-angle distribution of rj mesons with respect to the reaction plane for different 
bins in multipUcity Mreact (reaction centrahty) (a) - (c) for the experiment '^^Ni+'^^Ni at 1.9 AGeV 
and (d) - (e) for the experiment ^°Ca+'*''Ca at 2 AGeV, see Table |[ The corresponding values of 
the mean number of projectile-like spectator nucleons (Asp) are indicated. 
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Figure 23: The parameters Wj'^"'^ for ij mesons as a function of transverse momentum pt (av- 
eraged over multiplicity Mreact)'- boxes indicate results from the experiment ^^Ni+^^Ni at 
1.9 AGeV and triangles the results from the experiment '"'Ca+"'^*Ca at 2 AGeV. 

reconstructed reaction plane with respect to the true one, see sect. 3.2. The resulting 
dependence of 112''"'^ as a function of transverse momentum pt is shown in Fig. E3^ - p4 for 77 
and tt" mesons for both the Ni+Ni and Ca+Ca collisions. Note that the error-bars shown 
in Fig. E3^ - G4 represent statistical errors only. The systematical errors in the resulting 
parameter V2^"'^ are dominated by two sources: 

i) the uncertainty in the determination of the reaction-plane resolution {cos2/S.(j)pi) ( see 
sect. 3.2), which leads to a relative error in w|'''"^ of the order of 15%; 

ii) the uncertainty induced by the variation of the azimuthal anisotropy of the combina- 
torial background with ^2^ ~ 0.01 which was estimated from a Monte-Carlo simulation and 
leads to relative errors in ti|'""'^ of approximately 2%. 

Other sources of systematical errors like detector non-uniformities and occupancy of the 
TAPS spectrometer are found to be even smaller. 



5 Discussion 

The data shows strong negative elliptic flow of 77 mesons ( Wj'""^ < ) for all studied bins in 
transverse momentum pt and reaction centrality, indicating a preferred emission of 77 mesons 
perpendicular to the reaction plane. In contrast, the v*2^'^ values for iv^ are close to zero 
within the error bars except for peripheral Ni-fNi collisions {{Agp) = 51), where a weak 
elliptic-flow signal is observed, changing from positive to negative sign with increasing pion 
transverse momentum pt, see Fig. M 

In order to compare the present results with previous measurements of elliptic flow 
(l|, |l|, |l|, |l|, nil (see Table |), we plotted the wf"'^ values for neutral and charged 
pions near midrapidity as a function of the beam energy, see Fig. E5a. It can be seen that the 
pion elliptic flow undergoes a transition from out-of-plane to in-plane emission around beam 
energy of 2 AGeV, which can be compared with theoretical prediction |33] . This effect can 
qualitatively be understood by taking into account the collision dynamics. The spectator 
nucleons leave the interaction region after the so-called passage time ipass of the order 2R/7/3 
(where R is the nuclear radius and 7 is the Lorentz-contraction factor) , and the interaction 
with the spectator matter of particles produced later are not significant pa, M pa, |62| . 
Recent RQMD model calculations for the Au-f An system at a beam energy of 2 AGeV [ p4| 
show that the time-averaged values of V2 for both pions are positive (indicating a preferred 
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Figure 24: The parameters vl'""'^ for 7r° mesons as a function of transverse momentum pi for 
different bins in multiplicity Mreact- (a) - (c) for the experiment ^^Ni-|-^^Ni at 1.9 AGeV and 
(d) - (e) for the experiment 4"Ca-l-""*Ca at 2 AGeV. 
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Figure 25: (a) Beam energy dependence of the elliptic flow signal ^2''"'^ for neutral (circles) 
and charged (triangles) pions and rj mesons (squares) near midrapidity. Open symbols 
correspond to data deduced from published anisotropics and full symbols to our data (see 
text); (b)The parameter v^^"^ for 77 mesons as a function of the number of spectators (Asp). 
The dash-dotted line represent results of the pure geometrical absorption model (see text). 



28 



A+A 


E (AGeV) 


TT 


central! ty 


pt (MeV/c) 


^meas 


(cos2A0) 


^U-ue 


Rcf. 


Ar+Au 
Ar+Au 


0.095 
0.095 




semiperiph. 

ZPLF = 11-12 

peripher. 
ZpLP=14-17 


< 200 

< 200 


-0.12±0.05 
-0.23±0.09 


0.44^ 
O.44I 


-0.27±0.11 
-0.52±0.21 




y 
y 


Ar+Al 


0.095 


TT^' 


scmiccntr.. 


- 


-0.20±0.08 


0.49^ 


-0.42±0.16 




141 


Bi+Bi 


0.4 


TT^ 


MUL2 
MUL3 


360-600 
360-600 


-0.15±0.02 
-0.13±0.02 


0.35^ 
0.412 


-0.43±0.06 
-0.32±0.05 




18 




Bi+Bi 


0.7 


n^ 


MUL2 
MUL3 


360-600 
360-600 


-0.12±0.01 
-O.lOiO.Ol 


0.492 
0.512 


-0.24±0.02 
-0.20±0.02 




18 




Bi+Bi 


1.0 


n^ 


MUL2 
MUL3 


360-600 
360-600 


-0.08±0.01 
-0.09±0.01 


0.50 2 
0.552 


-0.16±0.02 
-0.16±0.02 




18 




Au+Au 


1.0 


tt" 


semicentr. 


400-600 


-0.19±0.07 


0.71^ 


-0.27±0.10 




16 




Au+Au 


11.5 


n^ 


peripher. 


400-500 


- 


- 


0.01-0.02 




21 
23 


- 


Pb+Pb 


158 


7r± 


scmiccntr. 


400-600 


- 


- 


0.02±0.01 




Pb+Pb 


158 


TT+ 


scmiccntr. 


50-800 


- 


- 


0.05±0.02 




22 





Table 7: Elliptic flow of pions for different beam energies. Comments: ^ - correction 
factors {cos2A(j)) were estimated from the simulations using the published uncertainty of 
the reaction plane determination. 2 _ correction factors taken from ref. |23] 



emission in the reaction plane), but pions emitted at the time scale less than tpass —13.5 
fm/c (see Table pi) exhibit negative elliptic flow. It is known that the freeze-out time of 
pions is strongly dependent on their transverse momentum and that most of the high pt 
pions freeze out early in the reaction ||3|, ^, ^. Therefore, the observed nearly isotropic 
emission of neutral pions, except the high pt pions in light colliding systems at 2 AGeV, can 
qualitatively be understood. 

On the other hand, the observed strong negative elliptic flow of 77 mesons indicates that a 
significant part of them freezes out while the spectators are still close to the participant zone. 
The magnitude of elliptic flow is comparable with that for pions observed in heavy colliding 
systems at energies 0.1-.4 A GeV Q - Q, see Fig. |25|a. However, the pt dependence is 
different, see Fig. ^ Within large error bars, the w^'^"'^ values for 77 mesons do not vary 
with transverse momentum in contrast to the ones of pions |Q - [|8| . In order to estimate 
the effect due to "shadowing" by spectators one can apply the formula 



= 0.5(1 - R)/{1 + R) , with B = exp{L/\), 



(15) 



where A is the mean free path for rj mesons in cold spectator matter (A,, ~ 1-2 fm [^) and 

1/3 
L = 2 ■ Asp fm is the mean thickness of spectator matter. The resulting dotted line in 

Fig. pSp describes the data for 77 mesons fairly well. 

The approximation by a shadowing scenario through cold nuclear matter for describing 
the strong negative elliptic flow of 77 mesons is also supported by different microscopic model 
calculations, which show that 77 mesons freeze out at earlier times as compared with pions 
|1(]| , ^, 65 . As an example. Fig. ^ shows the prediction of the BUU model ^^ for the time 
evolution of the A, A^*(1535), 7/ and tt yields in comparison to the density of the central 
region for central ^°Ca-|-'**'Ca collisions at 2 AGeV. This figure shows that the rj mesons 
decouple from baryons much earlier than the pions. Therefore, the shorter freeze-out time 
may explain the much stronger azimuthal anisotropy observed for 77 mesons in comparison 
with pions. 

In summary, we have studied simultaneously the azimuthal angular distributions of tt*^ 
and 77 mesons emitted at midrapidity in the two colliding systems ^*Ni-|-^^Ni at 1.9 AGeV 
and *°Ca-|-"°*Ca at 2 AGeV. We observed a strong negative elliptic flow of r/ mesons. The 
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System 


R(fm) 


Energy (A GeV) 


/3(c) 


7 


%ass{fm/c) 


Ni+Ni 


4.48 


2.0 


0.72 


1.44 


8.7 


Ca+Ca 


3.96 


2.0 


0.72 


1.44 


7.7 
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6.74 


2.0 


0.72 


1.44 


13.5 
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6.74 


1.0 


0.59 


1.24 


18.5 



Table 8: The nuclear radius ( R=1.6-A^/'^ fm ), velocity j3 (in units of the speed of hght), 
Lorencz-contraction factor 7 and passage time ipass for the systems Au+Au at 1-2 A GeV 
and Ni+Ni, Ca+Ca at 2 A GeV. All quantities are evaluated in the nucleon-nucleon center- 
of-mass system. 
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Figure 26: The prediction of BUU model for the time evolution of various quantities in a 
central colhsion of ^OCa+^OCa at 2 AGeV |6). 

elliptic flow of tt^ mesons is very weak in contrast to data obtained for heavy colliding 
systems at 1 AGeV. The data can qualitatively be explained by the final state interaction 
of mesons in the spectator matter taking into account the collision dynamics. 

This work was supported in part by the Granting Agency of the Czech Republic, by the 
Dutch Stichting FOM, the French IN2P3, the German BMBF, the Spanish CICYT, by GSl, 
and the European Union HCM-network contract ERBCHRXCT94066. 
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Figure 27: Mass spectrum of charged particles detected by TAPS for the ^*Ni+^*Ni 
system. Clearly separated peaks corresponding to charged pions, protons and deuterons 
are seen. 



6 Appendix A: Flow of charged baryons 

As TAPS allows the identification of charged particles by time-of-flight discrimination, pulse- 
shape discrimination and the charged-particle veto counters, we used this feature to verify 
our method of the reaction-plane reconstruction. As an example. Fig. |2^ shows the mass 
spectrum of charged particles detected by TAPS for the reaction ^^Ni-|-^^Ni at 1.9 AGeV. In 
this figure clearly separated peaks corresponding to charged pions, protons and deuterons 
are seen. Fig. Pq shows the experimental azimuthal-angle distributions of identified charged 
baryons (protons, deuterons and tritons) detected in TAPS with respect to the reaction 
plane for both studied reactions. The directed flow of charged baryons in the target-like 
rapidity region (the preferential emission in the reaction plane) is clearly observed. This 
verification demonstrates the quality of the reaction-plane determination for both colliding 
systems. 
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Figure 28: Azimuthal-angle distribution of charged baryons (protons, deuterons, tritons) at 
target-like rapidity (0.3 < yiab < 0.5). The left part contains the ^^Ni+^^Ni data at 1.9 AGeV 
and the right part the 4"Ca+""*Ca data at 2 AGeV. 
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